High-field magnetoresistance revealing scattering mechanisms in graphene 
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We show that the type of charge carrier scattering significantly affects the high-field magnetore- 
sistance of graphene nanoribbons. This effect has potential to be used in identifying the scattering 
mechanisms in graphene. The results also provide an explanation for the experimentally found, 
intriguing differences in the behavior of the magnetoresistance of graphene Hall bars placed on 
different substrates. Additionally, our simulations indicate that the peaks in the longitudinal re- 
sistance tend to become pinned to fractionally quantized values, as different transport modes have 
very different scattering properties. 

PACS numbers: 72.10.-d, 72.80.Vp, 73.23.-b 



Despite lots of both experimental and theoretical ef- 
fort, the issue of charge carrier scattering in graphene 
has not been completely settled [Il-Q. On SiC>2 sub- 
strates, scattering by long-ranged Coulomb impurities 
has been proposed to be a major factor contributing to 
the conductivity near the Dirac points [![. Specifically, 
the measured minimum conductivity in samples with a 
large width to height ratio is explained by such scattering 
[H, 0- Also short-ranged scattering is thought to occur 
in graphene, caused by resonant scatterers, like adsorbed 
hydrocarbons or hydrogen atoms [1] or, in the case of 
narrow graphene nanoribbons, edge roughness 0]. As 
charge carrier scattering significantly affects the perfor- 
mance of graphene-based devices, knowledge of the scat- 
tering mechanisms is essential to realize the full potential 
of the material. 

A recent study found that the behavior of the high-field 
magnetoresistance of graphene on a SrTi03 substrate is 
clearly different from that on a SiC>2 substrate Q ■ Specif- 
ically, the magnetoresistance peak corresponding to the 
zeroth Landau level (LL) decreases as a function of de- 
creasing temperature on SrTi03 , whereas previous stud- 
ies have shown that it increases with decreasing temper- 
ature on SiOa [IM! ■ Although the reason for this qual- 
itative difference has not yet been found, the results lead 
to strong argumentation about the dominant scattering 
mechanism in graphene [3. 14|. 

In this Letter, we simulate high-field magnetotransport 
in graphene nanoribbons (GNRs) and present results that 
agree qualitatively with the measurements. As our cal- 
culations will show, the magnetoresistance of the zeroth 
LL depends strongly on the type of the disorder, imply- 
ing that high-field magnetoresistance measurements may 
be used to probe the scattering mechanisms in graphene. 
We also find that the magnetoresistance peaks at the 
other LLs tend to become pinned to fractional values of 
the resistance quantum, as the localization lengths of dif- 
ferent transport modes differ significantly. 

Our results are based on an orthogonal tight-binding 
model, with hoppings taken up to third nearest neigh- 



bors [la ]. The model system is a two-terminal device 
with semi-infinite leads. Compared with a Hall bar, a 
two-terminal model is computationally less expensive and 
gives more fundamental information about the scatterers, 
as the explicit design of the voltage probes is not taken 
into account. It has to be noted that models that take 
electron-electron interactions into account predict spin- 
polarized states [lit ], while our computations apply to 
systems with degenerate spins. 

The transmission is computed using standard Green's 
function based methods and the conductance is acquired 
from the Landauer formula. We simulate charged impuri- 
ties with Gaussian potentials V = ^4exp[— (r — ro) 2 /2a 2 ], 
with the parameter A selected to be +0.1 eV for half of 
the impurities and —0.1 eV for the rest, in order to pre- 
vent doping. The amplitude of the potentials roughly 
corresponds to the measured amplitude of the potential 
variation on the surface of Si02 [17]. As the dielectric 
constant of SrTi03 increases from 200-300 to 5000 when 
the temperature is lowered from room temperature to 
liquid helium temperature [3, EH, we select the range 
a of the potentials as a variational parameter. Shorter 
ranged potentials represent charged impurities on SrTi03 
at low temperatures, and longer ranged those at higher 
temperatures or on Si02- 

We start by comparing the scatterer resistances R s of 
a 500 run long and 15 nm wide armchair-edged GNR 
(AGNR) and a zigzag-edged GNR (ZGNR) in left hand 
side panel of Fig. Q] We acquire R s from the total re- 
sistance i?tot by using the expression R s = R to t — ^/M, 
where the last term is the contact resistance with M be- 
ing the number of transport modes. All resistances are 
given in units of h/2e 2 . Each ribbon contains Gaussian 
impurity potentials at a concentration of 1.3 x 10 11 cm~ 2 , 
which is of the same order of magnitude as an estimate 
based on scanning tunneling microscopy of graphene on 
Si02 [13] • The range of the potentials is set to 15 times 
the C-C distance a, and the strength of the perpendicu- 
lar magnetic field is 200 T. With wider ribbons, weaker 
magnetic fields are required to produce similar results. 
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FIG. 1. (Color online) Scatterer resistance R 3 for GNRs with 
long-range Gaussian potentials at B — 200 T with respect to 
the location of the zeroth LL Eo (left), and the band struc- 
tures of pristine ribbons (right). The upper row corresponds 
to a 15 nm wide AGNR, and the lower row to an equally wide 
ZGNR. The dashed lines in the R 3 plots indicate the positions 
of the LLs in the pristine ribbons. Resistance values predicted 
by Eq. {TJ for the maximum values of the peaks and Eq. Q 
for the shoulders of the peaks of the AGNR are also shown. 



The band structures of pristine ribbons are shown in the 
right hand panel. If the spin degeneracy is not taken 
into account, the band structure of the AGNR contains 
a single valley with doubly degenerate LLs, whereas the 
ZGNR exhibits two valleys each with nondegenerate LLs. 
In AGNRs, the degeneracy of the LLs is broken near the 
edges of the ribbon [2l[ ■ It can be seen that R s is zero be- 
tween the LLs, as the physical separation of the opposite 
edge channels, corresponding to charge carriers moving 
in opposite directions, prevents backscattering. At the 
LLs the charge carrier density extends over the central 
part of the ribbon, and peaks in the scatterer resistance 
are formed. 

A study of the resistance peaks reveals fractional quan- 
tization 



22j, |23J , due to significantly different scattering 
properties of the transport modes. While modest dis- 
order causes the inner channels to become almost com- 
pletely localized, the outer channels remain almost com- 
pletely conducting. The total number of channels that 
correspond to charge carriers moving in a specific direc- 
tion (left or right) at the nth LL in a GNR is given by 
2\n\ + 1, and in the case of the ZGNR, the innermost 
channels in both valleys are localized. Thus, for LLs 
with \n\ > 1, we obtain R s = l/(2|n| - 1) - l/(2|n| + 1), 



which simplifies to 

R, = 



An 2 - 1 ' 



Inl > 0. 



(1) 



Additionally, a resistance peak of 2/3 is found just below 
the zeroth LL, corresponding to two of the three channels 
being localized. In the case of the AGNR, on the other 
hand, all channels at the zeroth LL are localized, result- 
ing in a very high peak in the resistance. Furthermore, 
when \n\ > 0, two kinds of behavior are observed. On the 
inside (closer to half filling) of each LL, the two inner- 
most channels are localized, yielding the same quantiza- 
tion expression as for the ZGNRs, i.e. Eq. ([1]). However, 
on the outside of the LLs (farther away from half filling) , 
only the innermost channels are localized, leading to the 
expression 
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Because of this, in an AGNR, the peaks in R s tend to 
exhibit steps or shoulders when \n\ > 0. 

A comparison with experimental data given in Ref . 24 1 
shows that the measured peaks in the longitudinal resis- 
tance of a multiterminal Hall bar of micrometer size differ 
roughly by 30 % from the values given by Eq. (JXJ) . Specif- 
ically, the measured resistance at the n = 1 to n = 3 LLs 
is lower than the prediction by Eq. (fTJ) , which would cor- 
respond to a situation where the inner channels are not 
yet fully localized. Additionally, the two-terminal model 
does not take into account effects caused by disordered 
contacts. 

We next study how the range of the potentials af- 
fects the heights of the peaks in R s . We define the 
energy-averaged typical scatterer resistance as (R s ,typ) = 
l/(exp(log(T to t))> - 1/(M) 0, where T is the transmis- 
sion function, the overlined quantity is ensemble averaged 
over random impurity configurations and the brackets 
stand for averaging over a finite energy window, AE, cor- 
responding to a finite bias voltage of V = AE/e. Here we 
have set AE to 1 meV. A statistical study of the peak 
heights in a 15 nm wide AGNR is presented in Fig. [2] 
where (a) shows how (-R s ,typ) evolves in a 15 nm wide 
and 100 nm long AGNR as the range of the impurity 
potentials is gradually changed, keeping the concentra- 
tion constant at 6.7 xlO 11 cm -2 . Each data point is 
based on an ensemble of 50 different realizations of im- 
purity locations and transmission averaging over 20 en- 
ergy points within the bias window. With short-ranged 
potentials, the central peak, corresponding to the n = 
LL, is shorter than the neighboring ones, but it increases 
drastically as a is made larger, as the charge carriers be- 
come localized. The peaks corresponding to other LLs 
increase only toward the quantized values, as only the 
inner modes become localized. The odd behavior of the 
n = peak may be attributed to the chiral symme- 
try of the graphene lattice. According to Ref. [26[, in 
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the presence of disorder that preserves the chiral symme- 
try, the |n| > LLs become broadened, whereas the ze- 
roth LL remains dclta-function-like. In our simulations, 
modest concentrations of short-ranged scatterers or rel- 
atively short-ranged Gaussian potentials, that cause a 
clear broadening of the resistance peaks at the \n\ > 
LLs, do not lead to a significant broadening the zeroth 
LL, if compared with the experimentally relevant energy 
scale set by a small bias voltage. This indicates that 
modest disorder affects the chiral symmetry of graphene 
only weakly, whereas stronger potential disorder, simu- 
lated by longer-ranged Gaussian impurities, leads to a 
broadening of all LLs. 

Fig. [U (b) shows the energy-averaged typical scatterer 
resistances at n = and n = 1 as functions of ribbon 
length, for a constant concentration (6.7 xlO 11 cm" 2 ) of 
impurity potentials with range a = 15a, as well as for 
a constant concentration (1.7 xlO 12 cm" 2 ) of monova- 
cancies. In the presence of the potentials, the resistance 
at n — grows exponentially with length due to strong 
localization, whereas the resistance at n = 1 approaches 
the fractionally quantized value of 2/3, given by Eq. ((T|). 
On the other hand, when the disorder is formed by va- 
cancies, the n = LL has clearly smaller resistance than 
the n = 1 LL. As shown in the linear plot of Fig. [2] (c), 
the n = peak grows linearly, whereas the n — 1 peak is 
already close to the fractionally quantized value, indicat- 
ing strong localization of the inner modes. Although not 
shown in the figure, strong enough disorder may cause 
scattering also in the outer channels, and thus an unpin- 
ning of the fractionally quantized resistance. 

The energy-averaged resistance corresponds to a lo- 
calization length £ [27j of roughly 60 nm at the zeroth 
peak for potentials with a = 15a. This can be turned 
into a scattering cross section E of 2.5 nm by using the 
expression £ = WLM/N£ [HHi], valid forsystems be- 
longing to the unitary transport ensemble [30j, where N 
is the number of impurity potentials and W and L the 
width and length of the ribbon, respectively. The scat- 
tering cross section is strictly valid only for a 15 nm wide 
AGNR, as the arbitrarily located potentials extend partly 
outside the ribbon edges. When |n| > 0, we approximate 
that only M' of the totally M modes exhibit scattering, 
and express the total resistance in the localized regime as 
R^j. = M'exp(-L/£')+M-M', where £' is the localiza- 
tion length corresponding to the M' modes only. Fig. [3] 
(c) shows R' = exp(L/£')/M' for the lowest LLs with 
fj = 15a, indicating a clearly localized behavior, with 
£' = 15 nm and £' = 20 nm. The inner modes at n = 1 
show a larger scattering cross section than the modes at 
n = 0, but the almost perfectly conducting outer modes 
limit the resistance to a fractionally quantized value. As 
the localization lengths of the inner modes are very short, 
one has to be careful when interpre ting resistivities given 
by the Ohmic formula p = RL/W [311 ] . 

If the Gaussian potentials are replaced by monovacan- 
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FIG. 2. (Color online) Statistical results showing the heights 
of peaks in the energy-averaged typical scatterer resistance, 
corresponding to the lowest Landau levels, in a magnetic field 
of 200 T. (a) Peak heights as a function of the range of the 
potentials in a system with a constant concentration of Gaus- 
sian impurities. The lines are to guide the eye. (b) Peak 
heights as a function of ribbon length for both Gaussian po- 
tentials and vacancies. The magnetic field is 200 T. (c) The 
results for the vacancies shown on a linear scale, highlighting 
the linear growth of the n — peak, (d) Growth of {R' Byt y P ) 
for monvacancies when n = 1. 



cies, £' for the n = 1 LL decreases to 8 nm, while the 
resistance peak at n — grows linearly as R s = AE/TT, 
yielding S = 0.065 nm. This indicates a qualitative dif- 
ference between long- and short-ranged impurities, as 
the latter scatter charge carriers only marginally at the 
zeroth LL. When the range of the potentials is short 
enough, the n = peak is much smaller than the n = ±1 
peaks, but once the range increases, the n = peak grows 
without limits, whereas the growth of other peaks is lim- 
ited by the presence of modes that exhibit only very lim- 
ited scattering. 

Finally, we compare the effects of impurity disorder 
in high magnetic field with those in zero magnetic field, 
using a test system that also has short-ranged disorder. 
Fig. [3] shows (i? s ,typ) computed from an ensemble of 100 
AGNRs of length 200 nm and width 15 nm, both with 
and without long-ranged impurity potentials, and with 
vacancies at a concentration of 6.7 xlO 11 cm" 2 . With 
no impurity potentials, the peak corresponding to the 
zeroth LL is negligible. If bulk vacancies are replaced by 
edge disorder, similar results are acquired (not shown). 
When four impurity potentials (two positive and two neg- 
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FIG. 3. (Color online) Typical scatterer resistance -Rs,t yP of a 
120-AGNR with vacancies but no long-range disorder (dashed 
red/gray), and in the presence of both vacancies and long- 
ranged disorder formed by N = 4 Gaussian potentials (black) . 
The peak at the zeroth LL has been averaged over an energy 
window of 1 meV. The inset shows the corresponding average 
conductances with respect to the Dirac point Ed at B = T. 



two-dimensional graphene, the phenomenon of minimum 
conductivity will prevent the observation of localization 
effects 

To conclude, in graphene nanoribbons, the magnetore- 
sistance at the lowest Landau level depends strongly on 
the type of the disorder. At higher Landau levels, and at 
the zeroth Landau level in zigzag edged ribbons, the mag- 
netoresistance may become pinned to fractionally quan- 
tized values due to significantly different conductance 
properties of the inner and outer transport modes. In 
armchair-edged ribbons, the resistance peaks at higher 
Landau levels tend to exhibit shoulders. According to 
our results, a decrease in the range of impurity potentials 
will cause a significant decrease in the measured magne- 
toresistance at the zeroth Landau level, which is in line 
with recent experimental results on a SrTiC>3 substrate. 

This research has been supported by the Academy 
of Finland through its Centres of Excellence Program 
(project no. 251748). We thank Mikko Ervasti, Ilja 
Makkonen and Antti-Pekka Jauho for commenting on the 
manuscript. 



ative), with range a = 8a, are added to the system, a 
high resistance peak is formed at n = 0. This contrasts 
clearly against the behavior at zero magnetic field, shown 
in the inset. The ensemble averaged conductance of an 
AGNR at B = T shows that the addition of impurity 
potentials does not produce any significant change in the 
conductance of a GNR in the absence of a magnetic field. 

Our results are valid for phase coherent systems. When 
the ribbon is longer than the phase coherence length, lo- 
calization effects will increase with lowering temperature, 
until the ribbon becomes fully phase coherent. Unless al- 
ready pinned to quantized values, the heights of the peaks 
in R s are expected to increase when the temperature is 
lowered, if the system is longer than the corresponding 
localization lengths. Importantly, the n — peak of an 
AGNR does not become pinned, which may correspond 
to the behavior observed in Hall bar experiments on Si02 
T3~13j]. On the other hand, in the recent experiment on 
SrTi03 jH, a decrease of the temperature and thus an 
increase of the substrate dielectric constant produced a 
decreasing peak at the n — LL, whereas the heights of 
the peaks at the \n\ > LLs increased. According to our 
simulations, this is explained by a reduced range of the 
impurity potentials. At zero magnetic field, no change 
in conductivity was found with decreasing temperature 
0| , but our simulations indicate that this does not neces- 
sarily imply that charged impurities are absent. As the 
inset in Fig. [3] shows, at zero magnetic field scattering 
by potential disorder may be completely masked by scat- 
tering by short-ranged disorder. On the other hand, in 
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